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ABSTRACT 

In  the  paper,  we  consider  models  of  SIR  -  SI  dengue  epidemic.  The  model  considers  the  effect  of  loss  (partial)  of 
immunity  on  different  state  variables.  Therefore  the  effect  of  loss  of  immunity  or  the  immunity  parameter  (yt)  plays  a  very 
important  rule  in  dengue  epidemic  model  and  gives  the  possibility  of  occurrence  of  the  dengue  diseases. 
The  characteristic  roots  of  the  model  at  disease  -free  equilibrium  (DFE)  point  are  the  real  and  opposite  sign, 
which  indicate  that  no  occurrence  of  dengue  virus  infection  since  there  are  no  infected  human  or  infected  mosquito 
population.  At  DFE  point,  every  human  in  the  population  is  healthy  and  not  infected  with  the  virus.  The  characteristic 
roots  at  endemic  equilibrium  point  are  all  negative  (real  part )  and  complex;  indicate  that  the  focus  of  dengue  fever  would 
be  stable.  The  estimated  value  of  the  basic  reproduction  number  is  2.73  with  a  range  (2.14,  3.02).  The  source  for  the  data 
is  NVBDCP. 

KEYWORDS:  Dengue,  Endemic,  Immunity,  Basic  Reproduction  Number,  and  Stability 


Article  History 

Received:  28  Feb  2018  I  Revised:  12  Mar  2018  I  Accepted:  16  Mar  2018 


1.  INTRODUCTION 

1.1  Global  Scenario  of  Dengue 

Dengue  vims  infections  are  a  serious  cause  of  morbidity  and  mortality  in  most  tropical  and  subtropical  areas  of 
the  world  especially  in  Southeast  and  South  Asia,  Central  and  South  America,  and  the  Caribbean.  Dengue  is  the  most 
rapidly  spreading  vector-borne  viral  disease  in  the  world  (Figure  1).  Unlike  other  vector-borne  diseases,  it  is  transmitted 
from  infected  human  to  a  female  Aedes  Aegypti  mosquito  by  a  bite  and  is  the  main  vector  for  dengue.  The  World  Health 
Organization  (WHO)  in  1998  has  listed  dengue  as  the  tenth  leading  cause  of  death  among  all  infectious  diseases. 
Severe  dengue  (also  known  as  dengue  hemorrhagic  fever)  is  found  in  tropical  and  sub-tropical  locations  in  most  Asian  and 
Latin  American  countries.  An  estimated  500,000  people  with  severe  dengue  require  hospitalization  each  year,  a  large 
proportion  of  these  are  children.  About  2.5%  of  those  affected  die  (Kurane  and  Takasaki,  2001).  The  incidence  of  dengue 
has  grown  dramatically  around  the  world  in  recent  decades.  The  actual  numbers  of  dengue  cases  are  un  reported  and  many 
cases  are  misclassified.  One  recent  estimate  indicates  390  million  dengue  infections  per  year  (95%  credible  interval 
284-528  million),  of  which  96  million  (67-136  million)  manifest  clinically  (with  any  severity  of  disease) 
(Bhatt  et  ah,  2015).  Another  study,  of  the  prevalence  of  dengue,  estimates  that  3.9  billion  people,  in  128  countries,  are  at 
risk  of  infection  with  dengue  viruses  (Brady,  et  ah,  2012). 
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Member  States  in  3  WHO  regions  regularly  report  the  annual  number  of  cases.  The  number  of  cases  reported 
increased  from  2.2  million  in  2010  to  3.2  million  in  2015.  The  disease  is  now  endemic  in  more  than  100  countries  in  the 
WHO  regions  of  Africa,  the  Americas,  the  Eastern  Mediterranean,  South-East  Asia,  and  the  Western  Pacific.  In  2015, 
2.35  million  cases  of  dengue  were  reported  in  the  Americas  alone,  of  which  10,200  cases  were  diagnosed  as  severe  dengue 
causing  1,181  deaths  (WHO,  2016).  In  2013,  cases  have  occurred  in  Florida  (United  States  of  America)  and  Yunnan 
province  of  China.  Dengue  also  continues  to  affect  several  South  American  countries,  notably  Costa  Rica,  Honduras, 
and  Mexico.  In  Asia,  Singapore  has  reported  an  increase  in  cases  after  a  lapse  of  several  years  and  outbreaks  have  also 
been  reported  in  Laos.  In  2014,  trends  indicate  increases  in  the  number  of  cases  in  the  People's  Republic  of  China, 
the  Cook  Islands,  Fiji,  Malaysia,  and  Vanuatu,  with  Dengue  Type  3  (DEN  3)  affecting  the  Pacific  Island  countries  after  a 
lapse  of  over  10  years.  Dengue  was  also  reported  in  Japan  after  a  lapse  of  over  70  years.  The  year  2015  was  characterized 
by  large  dengue  outbreaks  worldwide,  with  the  Philippines  are  reporting  more  than  169,  000  cases  and  Malaysia  is 
exceeding  111,  000  suspected  cases,  giving  about  59.5%  and  16%  increase  respectively  as  compared  to  the  previous  year. 
Brazil  alone  reported  over  1.5  million  cases  in  2015,  approximately  3  times  higher  than  in  2014.  Also  in  2015,  India 
(Delhi  alone),  recorded  its  worst  outbreak  since  2006  with  over  15,000  cases.  (WHO,  2016) 


Figure  1:  Global  Distribution  of  Dengue/Dhf 

1.2  Scenario  of  Dengue  in  India 

The  scenario  of  dengue  outbreaks  in  India  has  been  recently  reviewed.  The  data  on  the  website  of  National  Vector 
Borne  Diseases  Control  (NVBDC)  program  show  that  dengue  has  been  endemic  in  16  states,  Andhra  Pradesh,  Goa, 
Gujarat,  Haryana,  Karnataka,  Kerala,  Madhya  Pradesh,  Maharashtra,  Punjab,  Rajasthan,  Tamil  Nadu,  Uttar  Pradesh, 
West  Bengal,  Chandigarh,  Delhi,  and  Puducherry  (Table  1).  During  2010-2012,  dengue  spread  to  the  remaining  states. 
Figure  2  shows  that  the  distribution  of  dengue  cases  among  the  states  of  India  by  2015  (see  also  Figure  3).  Although  the 
number  of  dengue  cases  has  shown  a  steady  rise  for  some  states  with  every  passing  year,  the  mortality  has  reduced. 
Compared  with  the  rest  of  South-East  Asia,  the  number  of  dengue  shock  syndrome  (DSS)  cases  in  India  remains  low 
(Cecilia,  2014).  Yang,  et  ah,  (2009a  and  2009b)  discussed  the  temperature  effect  on  the  Aedes  Aegypti  mosquito. 
Although  dengue  has  been  notifiable  in  India  since  1996,  the  disease’s  impact  has  been  underestimated  because  of 
insufficient  information  on  incidence  and  cost  of  dengue  illness.  Between  2006  and  2012  the  NVBDC  program  reported  an 
annual  average  of  20,474  dengue  cases  and  132  deaths  caused  by  dengue.  Regional  comparisons  suggest  that  these  official 
numbers  reflect  only  a  small  fraction  of  the  full  impact  of  the  disease  (Shepard  et  ah,  2014). 


Impact  Factor  (JCC):  4.1675 
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Feng,  et  al.,  (1997)  constructed  a  model  to  study  both  the  epidemiological  trends  of  the  disease  and  conditions  that 
permit  coexistence  in  competing  strains.  Table  2,  gives  dengue  outbreaks  in  India  (Delhi  was  frequently  affected  part). 


Table  1:  Distribution  of  Dengue  Cases  in  India  for  Some  Selected  States  during  2010-15 


Affected  States 

2010 

2011 

2012 

2013 

2014 

2015 

Andhra  Pradesh 

776 

1209 

2299 

910 

1262 

3159 

Assam 

237 

0 

1058 

4526 

85 

1076 

Bihar 

510 

21 

872 

1246 

297 

1771 

Goa 

242 

26 

39 

198 

168 

293 

Gujarat 

2568 

1693 

3067 

6272 

2320 

5590 

Haryana 

866 

267 

768 

1784 

214 

9921 

Karnataka 

2285 

405 

3924 

6408 

3358 

5077 

Kerala 

2597 

1304 

4172 

7938 

2575 

4075 

Madhya  Pradesh 

175 

50 

239 

1255 

2131 

2108 

Maharashtra 

1489 

1138 

2931 

5610 

8573 

4936 

Orissa 

29 

1816 

2255 

7132 

6433 

2450 

Punjab 

4012 

3921 

770 

4117 

472 

14128 

Rajasthan 

1823 

1072 

1295 

4413 

1243 

4043 

Tamil  Nadu 

2051 

2501 

12826 

6122 

2804 

4535 

Uttar  Pradesh 

960 

155 

342 

1414 

200 

2892 

Uttrakhand 

178 

454 

110 

54 

106 

1655 

West  Bengal 

805 

510 

6456 

5920 

3934 

8516 

Chandigarh 

221 

73 

351 

107 

13 

966 

Delhi 

6259 

1131 

2093 

5574 

995 

15867 

Puduchery 

96 

463 

3506 

2215 

1322 

771 

India 

28179 

18209 

49373 

73215 

38505 

93829 

Source:  NVBDCP,  June  2016 
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(Source  :http://www.nvbdcp.gov.in) 

Figure  2:  Distribution  of  Dengue/Dhf  India  Figure  3:  Situation  of  Dengue  from  1996-2010 


Due  to  the  lack  of  effective  drugs  and  vaccines,  dengue  is  considered  as  a  public  health  problem  that  carries  a 
huge  financial  burden  on  the  governments  also.  Therefore,  due  to  the  lack  of  effective  programs  to  control  the  disease, 
a  wide  range  of  dengue  models  was  developed  to  describe  and  characterize  the  dynamics  of  its  transmission 
(Isea  and  Puerta,  2014;  Rodrigues  et  al.,  2011;  Side  and  Noorani,  2013;  Nuraini,  et  al.,  2007  and  Nishiura,  2006). All  these 
authors  were  made  the  effort  of  proving  a  better  understanding  of  the  nature  and  dynamics  of  Dengue  disease  transmission. 
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Chowell,  et  al.,  (2007)  estimated  the  reproduction  number  from  spatial  epidemic  data  at  the  level  of  municipalities 
using  two  different  approaches,  (i)  using  a  standard  dengue  epidemic  model  and  assuming  pure  exponential  initial 
epidemic  growth  and  (ii)  fitting  a  more  realistic  epidemic  model  for  the  initial  phase  of  the  dengue  epidemic  curve. 

Various  aspects  of  dengue  infections  in  the  context  of  India  e.g.,  occurrence,  clinical  profile,  viral  isolations, 
serological  surveys,  pathogenicity,  and  vector  ecology  have  been  discussed  by  Pandya  (1982).  Modeling  of 
epidemiological  diseases  is  an  important  tool  for  understanding  disease  behavior  and  prosing  effective  strategies  in  fighting 
against  disease  spread.  Esteva  and  Vargas  (1998)  discussed  a  model  for  the  transmission  of  dengue  fever  in  a  constant 
human  population  and  variable  vector  population.  The  control  measures  of  the  vector  population  are  discussed  in  terms  of 
the  threshold  condition,  which  governs  the  existence  and  stability  of  the  endemic  equilibrium. 

In  the  present  paper,  we  consider  the  development  of  dengue  epidemic  model  with  loss  of  immunity  because  the 
recovery  from  infection  by  one  virus  provides  lifelong  immunity  against  that  virus,  but  confers  only  partial  and  transient 
protection  against  subsequent  infection  by  the  other  three  viruses.  There  are  4  closely  related  serotypes  of  the  virus  that 
cause  dengue  and  the  lifelong  immunity  developed  after  infection.  With  the  four  closely  related  viruses  that  can  cause  the 
disease,  there  is  a  need  for  a  vaccine  that  would  immunize  against  all  four  types  of  dengue  to  be  effective. 

2.  DEVELOPMENT  OF  MODEL 

The  dengue  epidemic  model  identifies  two  populations,  a  human  (host)  population  (Nh)  and  a  mosquito  (vector) 
population  (Nv).  The  human  population  has  three  state  variables;  the  people  who  may  get  infected  with  dengue  virus 
(Susceptible,  Sh),  people  who  are  infected  with  dengue  (infected,  Ih)  and  people  who  have  recovered  from  the  disease 
(removed,  Rh).  The  vector  (mosquito)  population  (Nv)  is  divided  into  two  compartments;  mosquitoes  that  are  potentially 
may  get  infected  with  dengue  virus  (Susceptible,  Sv)  and  mosquitoes  that  are  infected  with  dengue  virus  (infected,  Iv). 


Table  2:  Epidemiological  Studies  where  Dengue  Virus  was  Identified 


Year 

Region  where  Study  was  Conducted 

Type  of  Dengue  Virus  Detected 

1964 

Vellore,  Tamil  Nadu 

DV-2 

1966 

Vellore,  Tamil  Nadu 

DV-3 

1968 

Vellore,  Tamil  Nadu 

DV-  1,2,3  &  4 

1968 

Kanpur,  Uttar  Pradesh 

DV-4 

1969 

Kanpur,  Uttar  Pradesh 

DV-4  and  DV-2 

1970 

Hardoi,  Uttar  Pradesh 

DV-2 

1983 

Kolkata,  West  Bengal 

DV-3 

1985 

Jalore  town,  South-West  Rajasthan 

DV-3 

1988 

Delhi 

DV-2 

1990 

Calcutta,  West  Bengal 

DV-3 

1988 

Rural  and  urban  areas  of  Gujarat 

DV-2 

1993 

Mangalore,  Karnataka 

DV-2 

1996 

Ludhiana,  Punjab 

DV-  1,2,3  &  4 

1996 

Lucknow 

DV-2 

1996 

Delhi 

DV-2 

1997 

Delhi 

DV-1 

1996 

Delhi 

DV-2  (Genotype  IV) 

1997 

Delhi 

DV-1 

1996 

Rural  areas  of  Haryana 

DV-2 

2001 

Dharmapuri  district,  Tamil  Nadu 

DV-2 

2001 

Gwalior,  Madhya  Pradesh 

DV-2 

2001 

Chennai,  Tamil  Nadu 

DV-3 

2003 

Northern  India  (Delhi  &  Gwalior) 

DV-3 
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Table  2:  Contd., 


2005 

Kolkata,  West  Bengal 

DV-3 

2003 

Kanyakumari  district,  Tamil  Nadu 

DV-3 

2003-04 

Delhi 

DV-3  (subtype  III) 

2003-05 

Delhi 

2003  -  DV  -  1,2,3  &  4  2005  -  D  -  3 

2006 

Delhi 

DV-3 

2006 

Delhi 

DV-1  &  3 

2001-07 

North  India  (Delhi  and  Gwalior  region) 

DV-1  (Genotype  III) 

2006 

Delhi 

DV-1, 3  &  4 

2008 

Delhi  region 

DV-1, 2  &  3 

1956-2005 

Entire  country 

DV-2 

2002-06 

Delhi 

DV-1,  2,  3  &4 

2003 

Delhi 

DV-3  (Genotype  III) 

2008 

Ernakulam,  Kerala 

DV-2  &  3 

2007 

Rural  areas  of  Madurai,  Tamil  Nadu 

DV-3  (Genotype  III) 

2007 

Andhra  Pradesh 

DV-1  &  4  (Genotype  I) 

2003-08 

Different  parts  of  the  country 

DV-3  (Genotype  III) 

2007-09 

Delhi 

DV  1,  2,  3  &  4 

2009-10 

Pune,  Maharashtra 

DV-4  (Genotype  I) 

Source:  Gupta,  et  al.,  (2012) 


Therefore  the  model  is  SIR  -  SI  dengue  epidemic  model.  The  transfer  diagram  of  the  model  is  as  follows 


Human  Population 


Vi 
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'& 
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1  fa  1 

Figure  4:  Transfer  Diagram  of  the  Different  State  Variables  of  Dengue  Epidemic 


Where  ph  the  is  transmission  probability  from  Iv  (per  bite),  pv  the  transmission  probability  from  7h  (per  bite),  9i  the 
proportion  of  infected  individuals  who  arc  infectious,  l/ph  the  average  lifespan  of  humans  (in  days),  l/pv  the  average 
lifespan  of  adult  mosquitoes  (in  days),  yh  the  rate  of  recovery  of  infected  individual  and  Yl  the  rate  at  which  recovered 
individuals  lose  their  immunity. 


The  system  of  differential  equations  for  the  host  population  is 


dS 


dt 

dLh 

dt 


h  N  h  ~  h  ^IvSh  -  hSh+Yltf, 


h  1  l vS  h 
N  „ 


-  (  h+yj  i, 


dR  „ 

dt 


y a i h  -  (  h  +  Yi)  ih 


(1) 
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And  the  system  of  differential  equations  for  the  vector  population  is 


dS 


dt 


I  S 

M _ X _ 1  h  v  _  S 

'  x  1  V  -V  T  v  W  ’ 

N  „ 


dl 


iIhSv 


dt 


N 


(2) 


With  initial  conditions,  Sh  +  Ih  +  Rh  =Nh  =>  Rh=  Nh  -  Sh  -  Ih  and 


A  A 

Sv  t-  Iv - =>  Sv  -  Nv  -  Tv - Iv 

My  My 


also 


=[(Sh  ,Ih,RhSv,  Iv )  G  :  (Sn  +,  Ih  +  Rh  <Nh,Sv  +IV  < 

Using  (3),  the  system  of  differential  equations  for  host  and  vector  populations  can  be  rewritten  as 

dS, 


(3) 


dt 

dl 


h  =  N 

h1N  h 


T  S 

h  1 ivJh 


N, 


A  +ylRh 


h  _  h 


i  ivs, 


dt 

dR  , 
dt 

dl 

dt 


N 


(  h  +T/,)I  / 


(4) 


=  yJh  -(  h  +  y i  )i i 


,ihs 


N 


Some  assumptions  in  this  model  are  (i)  the  total  human  population  (Ah)  is  constant,  (ii)  there  is  no  immigration  of 
infected  individuals  into  the  human  population,  (iii)  the  population  is  homogeneous,  which  means  that  every  individual  of 
a  compartment  is  homogenously  mixed  with  the  other  individuals,  (iv)  the  coefficient  of  transmission  of  the  disease  is 
fixed  and  does  not  vary  seasonally,  (v)  both  human  and  mosquitoes  are  assumed  to  be  born  susceptible;  there  is  no  natural 
protection  and  (vi)  for  the  mosquito  there  is  no  resistant  phase,  due  to  its  short  lifetime. 


To  simplify,  we  normalize  the  system  of  equations  (4)  by  defining  new  variables 


P  = 


Rh  *  _  h 
—  and  s  =  — 
Nh  Nv 


b 

Al>b 


(5) 


Therefore,  using  (5),  the  system  of  equations  (4)  can  written  as 


V  =  P/d  1  -  POD)  -  «  P  (t)  s(t)  +  j  r  (t)  (6) 

q'  =  a  p  (t)s  (t)  -  (^  q(t)  (7) 

r  =  X  q(t)  —  )?2r(t)  (8) 

s'  =  y(l  -  -  Ss(t )  (9) 
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Where  after  reparameterization,  we  have 

Pi  =  Yh  +  fth.  P2  =  Yi  +  fth.  «  =  y  =  Mu  8  =  i*v,A  =  y„  and  ^  =  y. 
Solving  the  system  of  equations  (6)  -  (9)  for  the  state  variables  p,  q,  r  and  s,  we  get 

*  1  2  TM-/I  +  (  1  2  ~  ^1  ) 

(  1  2  ~  ^1  ) +  YM-/I  1  2 

P  TM-fcCl  — P  ) 


i(7M1-/)  +  (  i 


XX^ 


)  ) 


(10) 


r  = 


X  P*yM1  -p*) 


t  XX, 

1  2(7M1-P )+( 1 — L) ) 

2 


s  = 


TP„(1  -  p‘) 


*  AA , 

TP„(1-P  )  +  (  i - L)  , 

2 


The  system  of  equations  (6)  -  (9)  yields  two  equilibrium  points,  one  is  disease-  free  equilibrium  (DFE)  point 
P°  =  (1,  0,  0,  0)  and  the  other  is  endemic  equilibrium  point  P  =  (p  ,  q  ,  r  ,  s  ).  Where  p  ,  q  ,  r  and  s  are  defined  in  (10). 

Linearization  of  a  systems  of  differential  equations  (6)  -  (9)  on  the  DFE  point  (1, 0,0,0)  yields  the  following 
characteristic  equation 


f 


Vs  J 


0 

0 

0 


0 

X 

y 


A, 

0 

0 


0 

-5 


a  yp^ 

q 
o 


vsy 


Therefore 


-  b~X 
0 
0 
0 


0 

X 

y 


0 

-  2-x* 

0 


0 

-  -X* 


=  0 


(II) 


This  is  the  determinant  of  jacobian  matrix  at  DFE  point,  Pu  =  (1,  0,  0,0).  Solving  (11),  we  get  :  =- 
and  another  three  roots  are  obtained  by  solving  the  Characteristic  polynomial  equation. 
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+  (l+  2+)  +*12+1  +  2  +Ct  )  +(PlP2  +0-P2  )— 0  (12) 

Similarly,  linearization  of  the  system  of  equations  (6)  -  (9)  on  the  endemic  equilibrium  point,  P*  =  (p  ,  q  ,  r  ,  s*) 
gives  the  following  Jacobian  matrix  evaluated  at  P  using  (10) 


J(P  * )  = 


Sj_ 

11  s2 

s2 

0 

0 


Y(1  “  sT* 


B_ 


B_ 

~A 

0 

Q±_ 

Q2 


Where  Sj  —  (A  —  B) 


(13) 


A,  A, 

S2  =  TM- /,  (A  -  B)  +  A(  , - i-  ) 

2 

Q  1  =  W- 1,  (A  -  B) 


)) 

Q  2  =  A  J  Y|I  h  (A  -  B)  +  A(  , - L  ) 

2 

A  =  12  (7M-  h  1)  —  ^  i 

B  =  t  2  (7B  k  +  )  -  M,  1 

3.  EFFECT  OF  LOSS  OF  IMMUNITY 


We  consider  the  effect  of  loss  of  immunity  on  different  state  variables.  Here  the  recovery  of  an  individual  from 
infection  by  one  virus  provides  lifelong  immunity  against  that  virus,  but  confers  only  partial  and  transient  protection 
against  subsequent  infection  by  the  other  three  forms  of  dengue  viruses. 


jTime  in  ctays 


Figure  5:  Effect  of  Recovery  Rate  on  State  Variables,  when  Individuals  Confer  Permanent  Immunity 
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Figure  6:  Effect  Recovery  Rate  on  the  State  Variables,  When  Individuals  Lose  their  Immunity 

The  parameter  values  used  for  this  are  given  in  Table  3.  From  Figure  5,  it  is  observed  that  when  an  individual 
doesn’t  lose  his/her  immunity  the  maximum  number  of  individuals  can  be  protected  and  from  Figure  6,  it  is  observed  that, 
up  to  certain  time,  the  recovered  individuals  increase,  once  they  start  losing  their  immunity,  the  recovered  individuals 
move  to  susceptible  class  thereby  recovered  individuals  get  reduced. 

4.  NUMERICAL  ILLUSTRATION  OL  THE  MODEL 

The  following  table  gives  model  parameters,  their  values  and  the  sources  of  parameters. 


Table  3:  Parameter  Values  and  their  Sources 


Parameter 

Value 

Source 

Nh 

480000 

Rodrigues,  et  al.,  (2011):  Chowell,  et  al.,  (2007) 

0. 

4.5 

Estimated 

Pvh 

0.3 

Asmaidi,  et  al,.(2014):  Rodrigues,  et  al.,  (2011) 

Phv 

0.3 

Asmaidi,  et  al,.(2014):  Rodrigues,  et  al.,  (2011) 

Ph 

1/71  x  365 

Rodrigues,  et  al.,  (2011) 

'•,+T 

1/3 

Asmaidi,  et  al,.(2014) 

V* 

Rodrigues,  et  al.,  (2011) 

Fv 

1/11 

Rodrigues,  et  al.,  (2011):  Chowell,  et  al.,  (2007) 

Nv 

3  Nh 

Rodrigues,  et  al.,  (2011) 

To  determine  the  critical  point,  the  equations  (6)-(9)  were  set  equal  to  zero  and  using  the  parameters  value,  we  get 
0.00004(1  -  p)  -  4.5ps  +  0.25  r  =  0 

4.5  ps  -  0.33334  q  =  0  (14) 

0.33334  q  —  0.25  r  =  0 
0.15(1  —  s)q  -  0.1s  =  0 

The  characteristic  roots  at  D  F  E  point,  P°  =  (1,  0,  0,  0)  are  obtained  as 
A |  =*  -0.00004,  2  =  -0.25,  3  -0.9663  and  4  =  0.7330 

Similarly  the  characteristic  roots  at  the  endemic  equilibrium  point,  P  =  (p  ,  q  ,  r  ,  s  )  are  obtained  as 
A  ,  =  -0.3907,  A 2  ==-0.2741  +  0.1721/,  A.,  ==-0.2741  +0.1721/  and  A 4=  -0.0977. 

The  characteristic  roots  at  DFE  point  are  a  real  and  opposite  sign,  which  indicate  that  no  occurrence  of  dengue 
virus  infection  since  there  are  no  infected  human  population  or  infected  mosquito  population.  At  DFE  point,  every  human 
in  the  population  is  healthy  and  not  infected  with  the  virus.  The  characteristic  roots  at  an  endemic  equilibrium  point  are  all 
negative  (real  part)  and  complex;  indicate  that  the  focus  of  dengue  fever  would  be  stable. 
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5.  BASIC  REPRODUCTION  NUMBER 


An  important  measure  of  disease  transmissibility  is  the  epidemiological  concept  of  the  basic  reproduction  number. 
It  provides  an  invasion  criterion  for  the  initial  spread  of  the  virus  in  a  susceptible  population.  The  basic  reproduction 
number,  denoted  by  R0,  is  defined  as  the  average  number  of  secondary  infections  that  occurs  when  one  infective  is 
introduced  into  a  completely  susceptible  population.  Using  the  next  generation  operator  approach  (Diekmann  and 
Heesterbeek,  2000),  we  compute  the  basic  reproductive  number  R0  associated  with  the  disease-free  equilibrium  point. 
The  basic  reproduction  number  of  the  model  (4)  is 


Ro 


QjPhPv 

tnTwt+r/i) 


(15) 


We  observe  that  R0,  the  reproduction  number  of  dengue,  depends  on  the  mosquito  and  human  vital  parameters, 
on  the  fraction  between  the  susceptible  mosquito  and  the  total  human  population  size,  and  also  on  the  product  of  the 
transmission  coefficients  and  the  square  of  mosquitoes  biting  rate,  9f(lh[iv,  emphasizing  that  a  new  case  of  dengue  can 
occur  only  after  two  bites  from  the  same  mosquito.  On  the  other  hand,  we  also  observe  that  dengue  control  appears 
explicitly  in  the  expression  of  R0,  only  on  yh.  The  basic  reproduction  number  (  R0)  decreases  as  yh  increases  resulting  in  a 
decrease  of  mosquito  susceptible  population. 


An  equilibrium  point  P  is  biologically  meaningful  if  and  only  if  P  £  T,  where  T  is  defined  in  (3).  The  biologically 
meaningful  equilibrium  points  are  said  to  be  disease-  free  or  endemic,  depending  on  Ih  and  Iv:  if  there  is  no  disease  for 
both  populations  of  humans  and  mosquitoes,  that  is,  if  Ih  ~  lv  =0,  then  the  equilibrium  point  is  said  to  be  a  Disease-  Free 
Equilibrium  (DFE),  P°;  otherwise,  if  Ih  =£0  or  lv  =£  0,  in  other  words,  if  lh  >  0  or  Iv  >  0,  then  the  equilibrium  point  is 
called  endemic,  P* 


If  R0  <  1,  then  the  disease  cannot  invade  the  population  and  the  infection  will  die  out  over  a  period  of  time. 
The  amount  of  time  this  will  take  generally  depends  on  how  small  the  basic  number  is.  If  P0  >  1,  then  an  invasion  is 
possible  and  infection  can  spread  through  the  population.  Generally,  the  larger  the  value  of  R0 ,  the  more  severe, 
the  epidemic  will  be.  In  determining  how  best  to  reduce  human  mortality  and  morbidity  due  to  dengue,  it  is  necessary  to 
know  the  relative  importance  of  the  different  factors  responsible  for  its  transmission.  Figure  7  gives  the  uncertainty  of  the 
basic  reproduction  number  related  to  the  parameters  in  the  model. 

The  basic  reproduction  number  is  most  often  estimated  from  data  in  the  early  epidemic  phase,  that  is,  prior  to  the 
introduction  of  initial  interventions,  including  behavioral  changes,  and  at  a  time  when  the  effects  of  susceptible  depletion 
are  negligible.  Using  approximation  (Favier  et  ah,  2006)  to  (15),  we  derive  the  basic  reproduction  number  with  force  of 
infection.  A,  which  is  obtained  assuming  the  cumulative  number  of  infections  follow  exponential  (that  is,  cumulative  total 
oc  eAt  )  as 

K0  =  (i+-4—)(1+-)  (16) 

V  Hv  +Yh'  Hv 

Equation  (16)  gives  the  relation  between  the  basic  reproductive  number,  R0  and  the  force  of  infection,  A. 

Substituting  these  parameter  values  (Table  3)  (no  vector  control  was  in  course  during  this  epidemic)  in  equation 
(16)  we  obtain  the  estimated  value  of  R0  =  2.73  (with  minimum  and  maximum  values  relative  to  the  human  parameter’s 
range  of  R0  =  2.14  and  3.02)  with  an  initial  growth  rate  of  the  epidemic  (or  the  force  of  infection),  0.099.  We  have  used 
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dengue  cases  in  India  (1996-2015)  for  estimating  R0.  The  data  we  used  for  estimation  is  taken  from  national  vector  -borne 
disease  control  (NVBDC)  programme. 


Basic  reproduction  number- 1 


Figure  7:  Basic  Reproduction  Number  with  Varying  Infection  and  Recovery  Rates 

Figure  7  indicates  that  there  is  no  advantage  of  having  with  varying  recovery  rate  when  infection  rate  is  also 
varying.  We  observe  that  the  curves  similar  when  infection  is  varying,  where  the  constant  rate  recovery  is  adopted  and 
having  a  varying  recovery  rate  when  a  rate  of  infection  is  constant. 

6.  CONCLUSIONS 

The  characteristic  roots  at  DFE  point  are  a  real  and  opposite  sign,  which  indicate  that  no  occurrence  of  dengue 
virus  infection  since  there  are  no  infected  human  or  mosquito  population.  At  DFE  point,  every  human  in  the  population  is 
healthy  and  not  infected  with  the  virus.  Whereas  the  characteristic  roots  at  an  endemic  equilibrium  point  are  all  negative 
(real  part)  and  complex  with  the  negative  real  part,  indicates  that  the  focus  of  dengue  fever  would  be  stable. 

The  model  considers  the  effect  of  loss  (partial)  of  immunity  on  different  state  variables.  Because  the  recovery  of 
an  individual  from  infection  by  one  virus  provides  lifelong  immunity  against  that  virus,  but  confers  only  partial  and 
transient  protection  against  subsequent  infection  by  the  other  three  forms  of  dengue  viruses.  Therefore  the  effect  of  loss  of 
immunity  or  the  immunity  parameter  (y, )  plays  a  very  important  rule  in  dengue  epidemic  model  and  gives  the  possibility 
of  occurrence  of  the  dengue  diseases.  Table  2  gives  occurrence  of  different  forms  of  dengue  in  India.  For  the  sensitivity  of 
the  parameter,  we  can  observe  Figure  5  and  6.  The  estimated  value  of  the  basic  reproduction  number  is  2.73  with  a  range 
(2.14,  3.02). 
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